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ABSTRACT 

M22 is an anomalous globular cluster that hosts two groups of stars with 
different metallicity and s-element abundance. The star-to-star light-element 
variations in both groups, with the presence of individual Na-0 and C-N anticor- 
relations, demonstrates that this Milky- Way satellite has experienced a complex 
star-formation history. We have analysed FLAMES/UVES spectra for seven stars 
covering a small color interval, on the reddest horizontal-branch (HB) portion of 
this cluster and investigated possible relations between the chemical composition 
of a star and its location along the HB. Our chemical abundance analysis takes 
into account effects introduced by deviations from the local-thermodynamic equi- 
librium (NLTE effects), that are significant for the measured spectral lines in the 
atmospheric parameters range spanned by our stars. We find that all the anal- 
ysed stars are barium-poor and sodium-poor, thus supporting the idea that the 
position of a star along the HB is strictly related to the chemical composition, 
and that the HB-morphology is infiuenced by the presence of different stellar 
populations. 
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1. Introduction 

The horizontal-branch (HB) morphology in globular clusters (GCs) is mainly governed 
by metallicity (e.g. Arp et al. 1952), which therefore is considered as the first parameter 
regulating the HB stars position along the color- magnitude diagram (CMD). The more metal- 
poor the GC, the more extended is the HB. The so called "second-parameter problem" 
recognises that the reality is more complex as some GCs with similar metallicity exhibit 
different HBs. Hence, at least one second parameter is required to properly caracterize the 
HB morphology in GGs (e.g. Sandage & Wallerstein 1960, van der Berg 1965). 

While the second parameter has been widely studied since the sixties, its identity is still 
unclear. Stellar mass, rotation, cluster central density, cluster mass, helium abundance, and 
age, among others, have been all suggested as possible solutions but none of them is fully 
satisfactory and possibly more than one parameter is at work (see Catelan 2009, Dotter et 
al. 2010 and references therein). 

It is now widely accepted that nearly all the normal GGs host stellar populations with 
different light-elements abundance (Garretta et al. 2009, and references therein). Observa- 
tional studies have shown that in GGs co-exist stars belonging to a first generation, whose 
O and Na composition is compatible with field stars at similar metallicities, and to second 
generation(s), e.g. those manifesting a chemical pattern (0-dcpletion and Na-cnhancement) 
compatible with a formation from material that underwent H-burning at high temperatures. 
Since the same material depleted in O and enhanced in Na is expected to be also He- 
enhanced, these findings provide new insight to the second-parameter problem, and seem to 
support the pioneering studies by Norris & Freeman (1982) suggesting that He abundances 
may take a role in the distribution of stars along the HB of GGs (e.g. D'Antona et al. 2005, 
Lee et al. 2005, Piotto et al. 2005, 2007, D'Antona & Galoi 2008). 

In this context, the GG M4 could be regarded as the Rosetta stone to connect multiple 
populations observed in the GMD with the HB morphology. This cluster indeed exhibits a 
strong dichotomy in the abundances of sodium and oxygen with the two main populations of 
Na-poor/O-rich first-generation, and Na-rich/O-poor second-generation stars defining two 
separate sequences along the red-giant branch (RGB) in the U versus U — B CMD. Due to 
the large differences in the CN and CH bands affecting the ultraviolet region of the spectrum, 
Na-rich stars define a sequence on the red side of the RGB, while Na-poor stars populate a 
bluer, more spread sequence (Marino et al. 2008). Interestingly, M4 hosts a bimodal HB, 
which is well populated both on the red and the blue side of the RR-Lyrae gap. It turns out 
that HB stars exhibit a bimodal Na and O distribution similar to that found along the RGB, 
with red-HB stars having solar [Na/Fe] and blue-HB ones being 0-depleted and Na-enhanced 
(Marino et al. 2011a, Villanova et al. 2012). The position of first-generation HB stars, whose 
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O and Na composition is compatible with field stars at similar metallicities, populate the 
red-HB. On the other hand, second-generation stars, supposed to form from material that 
underwent H-burning at high temperatures, distribute on the blue HB (see also Norris et al. 
1981 and Smith & Norris 1993). 

A similar scenario apphes for other GCs. Villanova et al. (2009) analysed spectra of 
seven HB stars with effective temperatures 8000< Teff <9000 K in NGC 6752 and found 
that six of them have a chemical content that resembles the field-halo composition, i.e. being 
the progeny of first stellar generation identified on the RGB. Successively, the connection 
between stellar populations and HB morphology in GCs has been confirmed for NGC 2808 
(Gratton et al. 2011), NGC 1851 (Gratton et al. 2012a), 47Tucanae (Milonc et al. 2012a, 
Gratton et al. 2012b), NGC 6397 (Lovisi et al. 2012), and M5 (Gratton et al. 2012b). 

While the majority of GCs are almost homogeneous in iron, recent studies reveal that 
there exists a group of anomalous GCs (AGCs), characterised by internal variations in the 
bulk heavy-element content, including iron and elements associated to slow neutron-capture 
(s-elements) . The chemical analogies with the two anomalous cases a;Centauri and M54, 
which are usually associated to the remnants of a dwarf galaxy cannibalised by the Milky 
Way, suggest the fascinating idea that AGCs could be the survived nuclei of more massive 
systems. 

M22 is certainly the prototype of an AGC. Indeed it hosts two groups of stars with 
different content of iron and s-elements, and this bimodality is associated to a different 
C+N+0 content (Marino et al. 2009, 2011b, 2012, Da Costa et al. 2009, Alves Brito et al. 
2012). Moreover, each individual group of s-rich and s-poor stars exhibits large star-to-star 
light-elements variation, with the presence of Na-0 and C-N anticorrelations (Marino et al. 
2011b). 

In this paper we analyse the chemical abundances for seven HB stars of M22 and 
investigate the possibility that stellar populations previously identified in this cluster could 
be related to a different position along the HB. 

2. Observations and data analysis 

Our dataset consists of UVES spectra collected within a project devoted to the study 
of sub-giant stars in M22 with GIRAFFE, (sec Marino et al. 2012). The fibers feeding the 
UVES spectrograph were centered on seven HB stars distributed along the more luminous 
portion of the blue HB of M22. The UVES spectrograph was used in the RED 580 setting 
providing spectra of R~45000. The spectra have a spectral coverage of ~2000 A with the 
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central wavelength at 5800 A. All our target stars were observed in the same FLAMES 
configuration in four different exposures of 46 minutes plus one exposure of 26 minutes, for 
a total observing time of 210 minutes. The single exposures were merged together, after 
they were reportad at the rest wavelength. The typical signal to noise ratio of the combined 
spectra is in the range 70 < S/N <100. Data were reduced using UVES pipelines (Ballester 
et al. 2000), including bias subtraction, flat-field correction, wavelength calibration and sky 
subtraction. 

We derived radial velocities (RVs) by cross-correlating the observed individual exposures 
with a spectrum simulated with atmospheric parameters close to our HB targets. The ob- 
served/template spectrum matches, after the heliocentric correction, yielded for our sample 
a mean RV of —144 ±2 km s~^ (rms=6 km s~^). This value is in reasonable agreement with 
the values in the literature (e.g. —143 ± 1 km s~^, rms=9 km s~^from Marino et al. 2012; 
—146.3 ± 0.2 km s~^, rms=7.8 km s~^, from the value reported in the Harris catalog). The 
rms/ a/ (A^ — 1) (where A^ is the number of available exposures) of our RVs measurements 
for each star is ~0.6 km s~^. This value has been taken as an estimate of the internal error 
associated with our RV values. The RVs derived from different exposures of the same star 
agree within our errors. Hence, as our spectra were collected in a time interval of ~3 months, 
the presence of binaries with period up to ~150-200 days can be excluded. 

The RVs of all the observed stars lie within ±3 a of the average RV (see Tab. [T|, that 
is the condition required for a star to be cluster member. However, in the case of M22, 
the RV-based criteria used to assess the cluster membership may be weak. Indeed M22 is 
projected toward the bulge, and the fleld is kinematically complex. Since the probability of 
flnding a very metal-poor fleld star with kinematics compatible with the cluster is extremely 
low, the metallicities obtained for our stars (listed in Tab. [l| see Sect. 2.1 ), all consistent with 



the M22 metallicities (e.g. Marino et al. 2009), constitute a more efficient tool to clean our 
sample from field contaminants. The status of He-burning star, as indicated by the derived 



atmospheric parameters (see Sect. 2.1) not compatible with a main sequence star (except 
maybe for one target), can be used as an additional proof, because the probability of a field 
contaminant with such properties is even lower. So, we conclude that all our observed stars 
are cluster members. 

We provide an estimate of the projected rotational velocities (wsinz) for our stars by 
deconvolving the full width at half maximum of the cross-correlation profile for the widening 
of the synthetic template computed with the instrumental profile of our data. Although 
a study of vsini is beyond the aims of this work, our values provide an estimate of the 



additional broadening to be applied to the lines when using spectral synthesis (see Sect. 2.1 ) 



The vsini of our sample stars ranges from 12 to 31 km s ^ in full agreement with values 
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determined for HB stars with similar T^s in a sample of other GCs (e.g. Recio-Blanco et al. 
2004; Behr et al. 2000). 

Our photometry included Johnson B, V, I bands from Peter Stetson's databas^ Pho- 
tometric catalogs were corrected for differential reddening as in Milone et al. (2012b) and 
were already used in Marino et al. (2011b), to which we refer the reader for details. The V 
versus B — I color-magnitude diagram (CMD) corrected for differential reddening is shown 
in Fig. [T| where we have marked the spectroscopic targets with cyan dots. 

Alpha-enhanced isochrones from the BaSTI and the Padova database (Pietrinferni et al. 
2004, Girardi et al. 2002) have been used to fit the observed CMD. We assumed reddening 
E(i? — l^)=0.35, [Fe/H]=— 1.85, and distance modulus (m — M)v=13.65, in agreement with 
values in the literature (Harris 1996, updated as in 2010). 

2.1. Atmospheric parameters 

Spectroscopic measurements were conducted with the local thermodynamic equilibrium 
(LTE) code MOOG (Sneden 1973) and by using model atmospheres interpolated from the 
grid of Castelli et al. (2004). In general, except for the cooler stars, too few Fe lines could be 
measured in our spectra to properly constrain all the atmospheric parameters. Therefore, 
we derived gravities (log g) from the isochrone that best fit our CMD, and associated to each 
star the log g value corresponding to the point on the isochrone with the smallest distance 
from the star. The distance is computed as in Gallart et al. (2003, Sect. 4) by enhancing 
the difference in color by a factor of seven with respect to the difference in magnitude. As 
pointed out by Gallard and collaborators, the factor seven, that is determined empirically, 
provides a weight to the color of a star that is larger than that of the magnitude. This takes 
into account the fact that, for a given isochrone, a star has a color that is better constrained 
than the magnitude, as any uncertainty in distance, gravity, or reddening, corresponds to a 
magnitude variation that is larger than that in color. 

Once log 5' had been fixed, effective temperatures (Teg) and microturbolence values (^t) 
were derived from the spectra, by satisfying the ionisation equilibrium between Fe I and Fe II 
abundances, and by removing trends in Fe abundance^ with reduced equivalent width (EW), 
respectively. We have used the non-LTE (NLTE) spectral code MULTI (Carlsson 1986; see 
Lind et al. 2012 and Bergemann et al. 2012, for a description of the NLTE modelling) to 



' http://www3.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/community/STETS0N/index.html 

^Fe I or Fe II, depending on the number of lines available for the different ionisation stages. 
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estimate the offset between Fe I and Fe II abundances expected for each star. It turns out 
that in LTE approximation the abundances from Fe I lines are lower than the ones from Fe II 
lines by ~0. 25-0. 28 dex, due mainly to large positive NLTE effects on Fel lines. This offset 
has been considered in the ionisation balance for the spectroscopic determination of Tcfj. 
Our approach has the important advantage that it is independent from photometry and is 
not affected by the differential reddening effects present across the face of M22. On the other 
hand, log g is more robustly internally fixed from isochrones, given that our stars distribute 
on an almost horizontal section of the HB and are expected to show similar gravities. For 
one star (#2562) no Fe I line could be measured from the spectrum, and we adopted the 
Teflf derived from the best-fitting isochrone. 

The final adopted atmospheric parameters for each star are listed in Tab. [T] Internal 
errors in log g and Tes depend on photometric uncertainties and the ionisation equilibrium, 
respectively. For our relatively bright stars, internal errors in the color {B — /) and in the 
magnitude V, including uncertainties due to the differential reddening correction, are 0.025 
and 0.015 mag respectively (Milone et al. 2012b). A shift of our target stars in the position 
along the ZAHB by these photometric uncertainties corresponds to a variation in log g of 
~0.02, that is so small to translate in negligible internal errors in the derived chemical 



abundances (see Sect. 2.2). To determine internal errors in temperatures, we adopted the 
following approach to all the stars in our sample: (i)we added the standard error of the mean 
for Fe I and Fe II in quadrature; (ii) then adjusted Teff until the quantity Fe I— Fe II, after 
the NLTE corrections, was equal to this value. The mean difference between the new T^s 
and the original values is ~170 K, that is an estimate of the internal uncertainties affecting 
our Teff determinations. For we measured the formal uncertainty in the slope between the 
Fe abundances and the reduced EWs, and adjusted until the formal slope was equal to this 
value. Given the relatively low number of spectral lines available, internal uncertainties in 
are large, and are on average equal to ~0.5 km s~^. 

For a comparison between HB and RGB chemical abundances, that is the subject of our 
study, an estimate of possible systematic errors affecting the model atmospheres is required. 
External errors are much harder to estimate. A comparison of atmospheric parameters de- 
rived with different techniques may be indicative of possible systematics affecting our adopted 
values. In Tab. [T| along with the adopted parameters we also show the Teg determined from 
isochrones, and resulting from the Tgff-^t empirical relation derived in Pace et al. (2006) 
from the least square fit of the HB data analysed by Behr et al. (2000). A comparison of the 
two sets of Teff suggests that the agreement is satisfactory for the three coldest stars (e.g. 
#2212, #2346, and #2548), whose difference among spectroscopic and photometric T^s is at 
most ~100 K. For the two hottest stars (#2297 and #90) spectroscopic Teff are significantly 
lower, suggesting that a possible systematic between the two sets of T^s values could be 
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higher for warmer stars. On average the difference between the two sets of T^s values (#166 
and #2562 have been excluded in this calculation) is Tcff^adoptcd— Teff4sochroncs=— 174 ± 90 K, 
that could be taken as a rough estimate of possible systematics affecting our Teg determi- 
nations. We note that, although the limited number of stars presented here do not allow us 
to provide a more exhaustive analysis of possible systematics, we consider the spectroscopic 
Tgff more reliable since the photometric ones are inconsistent with the Fe ionisation balance 
in NLTE. The difference between the adopted and the ones obtained by the Pace et al. 
relation is ^t.adoptod— ^t,Pace=— 0.28 ± 0.13 km s^^. To estimate the error provided by the 
adopted stellar models, we compared the values of Tefr and log g from the two independent 
set of isochrones used in this paper. Differences in temperatures and gravities are typically 
~80K and ~0.04 dex, respectively, with the BaSTI isochrones providing smaller values of 
Tgff and log g. In the following, we will use the log g provided by the Padova isochrones only. 

Our gravities values are based on the assumption that our stars lie on the ZAHB. 
During the He-burning phase, however, the stars are expected to decrease their gravity, 
which are about 0.2 dex lower for terminal-age HE stars. The photometric changes during 
this evolution should be negligible. Hence, log g values derived from isochrones could be 
overestimated by up to ~0.2 dex if the target is close to the end of its HB life. Moreover, 
canonical models with no He-enhancement, as those we used here, predict log 5' higher than 
logg from He-enriched models. In our models gravities derived from isochrones at Y=0.30 
are higher by ~0.25 dex than those derived by assuming primordial He. Some evidence for 
gravities systematically lower by 0.3 dex than those predicted by stellar models, have been 
shown in the case of u Gen (e.g. Moni Bidin et al. 2011, Moehler et al. 2011), but not in 
normal GCs (Moni Bidin et al. 2007, 2009). As M22 shares many chemical features with 
the most extreme GC u Gen (e.g. Marino et al. 2009; Da Gosta & Marino 2011), the current 
models could be similarly inadequate to describe HB stars in M22. For these reasons, gravity 
estimates independent on isochrones are important to have an idea of possible systematics 
in log g introduced by these effects. 

Two Balmer lines lie in the spectral range of our data: the H^ (~6563 A) and the H^ 
(~4861 A). As, in the Tgff range of our stars, the hydrogen lines are not very sensitive to the 
temperature, we can provide estimates of gravity from a minimisation of the observed line 
wings with theoretical spectra. To this aim we calculated spectral models for the adopted Tgff 
(Tab. [1]) using the spectral synthesis code SYNTHEQ Internal errors associated with these 
measurements are strongly dominated by normalisation uncertainties, that are much larger 
in the case of the H^ lying in the proximity of the blue-edges of the spectra. For this reason. 



- http://wwwuser.oat.ts.astro.it/castelli/sources/synthe.html 



- 8 - 



we preferred to use only the Hq, for the log g estimates. As an example, in Fig. [2] we represent 
the spectral synthesis for the star #90. Together with the model providing the minimum x^, 
at log (7=3. 13, we plotted the models at log (^=±0. 30 dex around the best-fitting model. The 
fit with synthetic spectra at different logs' suggests that, due to the quality of our observed 
spectra and the sensitivity of the Hq, line wings to log g, internal errors associated with these 
measurements are ~0. 15-0.20 dex. The logs' values derived from the H„ best-fit models, 
listed in Tab. [l| are all (apart for the star 7^2562) lower than those derived from isochrones. 
On average the difference among the gravities obtained from the two different techniques is 
—0.14 ± 0.08 dex (rms=0.17), that is an estimate of possible systematics affecting our logs 
determinations. We verified that lowering gravities by this quantity refiects in temperatures 
lower by 110 ± 8 K needed to establish the ionisation equilibrium. 

In the following, we will use log g provided by isochrones, because they provide a higher 
level of internal precision, and Tes from the ionisation equilibrium. For comparison purposes, 
atmospheric parameters obtained with different techniques are listed in Tab.[T} As discussed, 
our analysis suggests that the offsets with respect to the adopted parameters are at most of 
~200 K in Tefj, ~0.15 dex in logs, and ~0.30 km s~^ in The impact of these possible 



systematics on chemical abundances is discussed in Sect. 2.2 



Note that for one star in our sample (#166) the model atmosphere entirely derived from 
isochrones gives a large difference (0.68 dex) among Fe I and Fe II abundances. By deriving 
Teff and from the Fe lines, once log g has been fixed from isochrones (as done for the other 
stars), we need to significantly lower Teg and to 7450 K and 1.0 km s~^, respectively. 
While Teff obtained from the ionisation equilibrium is comparable with the lower-Teff stars 
in our sample, C,t is too low. Once fixed the other parameters, gravity from results slightly 
higher than the isochrone value (see Tab. [T|. The spectrum of this star shows a relatively 
sufficient number of Fe I (19) and Fe II (4) lines to attempt to entirely constrain parameters 
from spectral lines, e.g. to derive both Tog and log g we impose the excitation and ionisation 
equilibrium, respectively. The log g value that establishes the ionisation balance of iron 
lines, after the offset between Fe I and Fe II lines expected in the LTE approximation has 
been considered, is almost one dex higher than that provided by isochrone fitting, while we 
still get values lower than in the other stars. The reason for this discrepancy is unclear. 
Magnitude and color provided for star #166 by a different photometric catalog (from Monaco 
et al. 2004) are very close to those of this paper, thus suggesting that the high value of log g 
is physical. One possibility could be that #166 is an evolved Blue Straggler. Indeed, the 
comparison of the CMD of Fig. [T] with an 1.1 Gyr isochrone with the same metallicity as 
M22 shows that the color and magnitude of #166 matches the location of an evolved 1.8 
Mq star. We leave this issue to future investigations. Here, we use fully spectroscopically 
derived atmospheric parameters for this star, and suggest caution in considering its derived 
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chemical abundances. 



2.2. Chemical abundances 

Besides a few Fe I and Fe II absorption lines, our spectra show the Na resonance doublet 
at ~5890 A, and few lines of Mg, Ca, Ti, Ba, and Cr. Line broadening from stellar rotation 
is evident in the spectra, but even in the most extreme cases, the line profiles were close to 
Gaussian, so line EWs were measured by least-square fitting of Gaussian profiles to the data. 
Adopted log (7/ are from the NIST databas^for all the lines, except for Fe II for which we 
used the log gf from Melendez & Barbuy (2009). 

In the five (out seven) coolest stars in our sample we detect the two Ba lines at 
A ~4934 A, and ~6141 A. These two lines were measured by using spectral synthesis, to 
properly account for the blends with Fe lines affecting both lines and the effects of hyperfine- 
splitting. For the synthesis we used linelists based on Kurucz line compendium]^ apart from 
the Ba transition for which we added hyperfine structure and isotopic data from Gallagher 
et al. (2010). 

As Na resonance spectral lines are heavily affected by departures from LTE, we de- 
termined NLTE corrections to our abundances as in Lind et al. (2011), tailored with our 
Kurucz model atmospheres. The overall effect of the NLTE corrections to our abundances 
is a decrease of the derived [Na/Fe] abundances and of the line-by-line scatters. 

Our derived chemical abundances are listed in Tab. [2j For Fe and Na we report both 
the LTE and NLTE abundances. We note that the NLTE mean metallicity derived for 
HB stars, that is [Fe/H]= -1.65 ± 0.04 (rms=0.10) agrees within ~ 2a with the [Fe/H] 
derived on the RGB of the cluster [Fe/H]LTE = -1-76 ± 0.02 (rms=0.10) (Marino et al. 
2009). This agreement is satisfactory if we consider the different spectral lines measured in 
RGB and HB stars and the different analysis used in Marino et al. (2009, 2011b). Indeed, 
in the atmospheric parameters domain of the M22 RGB stars the NLTE corrections on Fe 
abundances are much smaller than for HB stars. As such, the ionisation equilibrium between 
Fe I and Fe II is not seriously affected by NLTE effects, and the atmospheric parameters 
were derived in LTE. Accounting for the small NLTE effects affecting Fe I lines in RGB stars, 
the mean metallicity of M22 would have been systematically higher by ~0.06 dex (Lind et 
al. 2012), thus decreasing further the HB-RGB mean metallicity offset. A more proper 



' http://physic5.ni5t.gov/PhysRefData/ASD/lines-form .html 
' http://wwwuser.0at.ts.a5tr0.it/castelli/linelists.html 
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comparison between RGB and HB stars in M22 has to take into account the different stellar 
groups hosted in the cluster, which is the subject of Sect. [3] 

To have an estimate of internal errors introduced by model atmospheres and EW uncer- 
tainties in the abundances, we repeated the abundance measurements by changing one at a 
time Teff, log g and by A(Teff)= ±170 K, A(log g)= ±0.02 dex, A(6)= ±0.50 km s-\ that 
are the internal uncertainties associated with the atmospheric parameters, as determined in 



Sect. 2.1 Metallicity has been changed by A([A/H])= ±0.10, as suggested by the derived 
dispersion in Fe II abundances, and EWs by A(EW)= ±8 mA, that is the typical error 
associated with our measurements as we have verified by comparing EWs obtained from 
single exposures of the same star. Internal uncertainties in chemical abundances have been 
determined by applying this procedure for all the stars in our sample. Then, the average un- 
certainties for each chemical species, in Tab. |3| have been taken as estimates of the errors in 
abundances. Assuming that the uncertainties listed in Tab. |3] are uncorrelated, we estimate 
total abundance uncertainties by summing in quadrature the various contributions. The 
resulting internal errors in chemical abundances (atotai in Tab. [3]) are typically of ~0. 10-0.20 
dex and, in most cases, are comparable to the observed dispersions. In Tab. [3] we also list the 
sensitivity of derived abundances to possible systematics that may affect our atmospheric 
parameters (see Sect. 2.1). 



3. Results 

The lower-left panel of Fig. [3] shows [Ba/Fe] as a function of Teg for the HB stars 
measured in this paper (circles) and RGB stars studied by Marino et al. (2009, 2011b, 
triangles). In these plots a correction of 0.06 dex in Fe has been applied to the RGB stars 



to take into account NLTE effects, as discussed in Sect. 2.2 



The lower panels of Fig. [3] show that for RGB stars, the barium distribution is clearly 
bimodal. We colored red and blue s-rich and s-poor stars, respectively, as they have been 
defined by Marino et al. (2009, 2011b). These color codes will be used consistently hereafter. 
The comparison of the histogram distribution of [Ba/Fe] for HB stars (lower-middle panel) 
and RGB stars (lower-right panel) reveals that all the analysed HB stars are consistent with 
being the progeny of s-poor stars. Despite that the errors associated with our HB [Ba/Fe] 
measurements are larger than the ones associated with RGB values, we note that: (i) the 
star showing the higher Ba abundance is the peculiar star #166, that we do not consider in 
the discussion; (ii) the mean HB Ba content of our analysed stars is low. Excluding the two 
stars #166 and #90 (for which we can provide only an upper limit to the Ba abundance), the 
mean [Ba/Fe] for HB stars is lower than the ones derived for s-poor stars on the RGB, but 
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still compatible with this value within a 3a level. Many factors can explain the discrepancy 
between the mean HB and s-poor RGB Ba abundance, such as the different lines analysed, 
NLTE effects, and/or systematics in our parameters scale. We note that the Ba line A4934 A, 
that was not used in the RGB analysis, gives Ba abundances systematically lower by ~0.15 
dex than the line A6141 A, that was instead measured on RGB stars. The sensitivity of 
[Ba/Fe] abundances to systematics in the model atmospheres listed in the three right-most 
columns of Tab. |3| suggests that, even accounting for the largest systematics, it is unlikely 
to reconcile the HB mean [Ba/Fe] with the s-rich RGB mean abundance. Therefore we can 
safely associate all the analysed HB stars to the s-poor stellar group. 

Additional arguments for the association between the s-poor group discovered on the 
RGB and the analysed HB stars come from the inspection of other chemical abundances. 
As discussed in Marino et al. (2011b), chemical properties characterising the two stellar 
groups with different s-process abundances include that (i) they have a different metallicity, 
with s-poor stars showing an overall metallicity ~0.15 dex lower than the s-rich stars, (ii) 
and neither the s-rich nor the s-poor group is consistent with a simple stellar population, as 
both groups exhibit significant variations in the abundance of light elements carbon, nitrogen, 
sodium, oxygen, and aluminium delineating individual C-N/Na-O/Na-Al (anti) correlations 
(see Marino et al. 2009, 2012). Unfortunately, the error on Fe for HB stars is much larger 
than for RGB stars. The mean [Fe/H] derived from the HB sample perfectly corresponds 
to the s-rich average metallicity, but it still agrees within a 3 a level with the mean value 
obtained for the s-poor giants. Therefore, we are not able here to use the relatively small 
difference in metallicity as a tool to associate the analysed HB stars to one or the other RGB 
stellar group. This is also because the possible interplay of systematics, such as the different 
analysed lines, different treatment of NLTE effects, can easily shift our mean metallicity. 

As sodium varies much more than [Fe/H] in the different M22 stellar populations, for 
our HB sample, [Na/Fe] abundances could provide some constraints on the belonging of these 
stars to the different stellar populations identified on the RGB. In the upper-left panel of 
Fig. |3| we plot [Na/Fe] against Tes for HB and RGB stars. Both s-rich and s-poor RGB stars 
span a large range in [Na/Fe] with s-rich stars having, on average, higher sodium content. 
The histogram distribution of [Na/Fe] for s-poor stars is multimodal, with about one half 
of the stars clustered around [Na/Fe]~ —0.20 dex (upper-right panel of Fig. |3]). None of the 
RGB s-rich stars have [Na/Fe]<0. The HB stars analysed in this paper span a much more 
narrow range in sodium. Their mean sodium content is [Na/Fe]=— 0.34±0.03 (rms=0.06) 
and they appear to belong to the same group of s-poor/Na-poor stars identified along the 
RGB. The mean Na abundance relative to Fe for these RGB stars is [Na/Fe]=— 0.13 ± 0.03 
(rms=0.09), which would decrease by ~0.06 if we account for NLTE corrections on Fe I RGB 
abundances. Despite the mean Na abundances of the analysed HB stars agrees within 3a with 
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the average value measured on s-poor/Na-poor RGB stars, we cannot exclude systematic 
effects due to the different lines analysed. In any case, the low mean [Na/Fe] measured for 
the HB analysed stars provides a further and more robust indication that the these stars 
belong to the s-poor group. In addition to this, within the s-poor group, they constitute the 
progeny of the Na-poor primordial population. 

Within observational errors the chemical abundances of the other analysed elements 
Mg, Ca, Ti, and Cr do not show evidence for internal variations. The a element abundances 
of Mg, Ca, and Ti relative to iron are typical of Population 11 stars, as also observed on 
the RGB (Marino et al. 2011), while chromium abundances relative to iron are roughly 
solar-scaled. 



4. Conclusions 

Observational evidence of a connection between the distribution of stars along the HB 
and their light element chemical content have been shown in the normal GCs M 4, 47 Tuc, 
NGC 2808 and NGC 6397 (Marino ct al. 2011, Villanova ct al. 2012, Gratton et al. 2011, 
Gratton et al. 2012, Lovisi et al. 2012). We have analysed spectra for seven stars in M22 
spanning a small range of temperature (7300< Tes <8300) along the HB, with the main 
purpose of measuring their Ba and Na abundances. Our results support the idea that the 
distribution of stars along the HB is somehow related to the multiple stellar population 
phenomenon also for the anomalous GC M 22. 

Recent studies on RGB and SGB stars have shown that M 22, differently from the most 
{normal) GGs, hosts two groups of stars with different content of iron, s-elements, and overall 
G+N+0 (Marino et al. 2009, 2011, 2012; Da Gosta et al. 2009, Alves Brito et al. 2012). 
Both the iron-/s-/CNO-rich and iron-/s-/CNO-poor group exhibits an internal scatter in 
the abundance of light elements (e.g. C, N, O, Na, Al) with the presence of Na-0 and 
G-N anticorrelations. Fe- and s-poor stars have, on average, lower sodium and nitrogen 
abundance than Fe- and s-rich ones. Indeed, in this cluster, the chemical tracing of multiple 
stellar populations along the HB must take into account not only light elements, but also 
s-process elements. All our analysed HB stars are consistent with having all the same Ba and 
Na abundances, suggesting that they are consistent with being the progenie of the Na-poor 
and s-poor first generation stars observed along the RGB. 

Similarly to what recently observed in normal GGs, the reddest HB stars of M 22 
belong to the primordial population, whose light-element chemical composition (e.g. Na) is 
consistent with field stars at similar metallicity. As a consequence of this, it is tempting to 
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speculate that He could be an important parameter in determining the HB morphology in 
GCs, as second generation(s) stars are supposed to born from H-burning processed material. 
He-rich stars should evolve faster than He-poor ones and hence, present-day Na-rich stars 
should be less massive than Na-poor stars. As a consequence of their mass, He/Na-rich and 
Hc-Na-poor stars should evolve into different HB regions, with Na-rich HB stars having also 
redder colors. However, direct and strong observational proofs supporting this idea are still 
missing. 

Some indication for a possible He enhancement of second gcncration(s) stars have been 
shown by Dupree et al. (2011) and Pasquini ct al. (2011), on the basis of the detection of a 
chromosphcric He line for Na-rich RGB stars. In addition to this, the comparison between 
the observed multiple sequences along the CMDs and theoretical models of a number of GCs 
(e.g. Piotto et al. 2007, Milone et al. 2012a) suggest that present-day Na-rich stars in GCs 
could be enhanced in He. On the other hand, the recent results of Moni Bidin et al. (2011) 
cast many doubts, as they show that blue HB stars in uj Centauri are not brighter than in 
other GCs, as would be expected as a consequence of a possible He enhancement. 

In M 22 a further complication comes from the presence of star-to-star variations in the 
overall C-l-N-l-0 content. Theoretical models, indeed, predict that the effective temperature 
of stars along the HB, as well as the mass loss along the RGB is affected by the CNO 
abundance (Cassisi et al. 2008, Ventura et al. 2009). The fact that all the five targets 
for which barium content has been measured show abundances consistent with the s-poor 
group identified on the RGB (Marino et al. 2009), showing C+N+0 abundances lower 
than in s-rich stars, supports the hypothesis that the CNO abundance can play a role in 
determining the distribution of stars along the HB. In general, the result that a small HB 
segment is populated by stars with almost the same chemical composition, supports the idea 
that the presence of multiple stellar generations with different chemical composition is a key 
parameter in determining the HB morphology. 

We warmly thank the anonymous referee for his/her comments, that helped to strengthen 
our results. We are grateful to Peter Stetson for providing photometry of M 22. APM ac- 
knowledges the financial support from the Australian Research Council through Discovery 
Project grant DP120100475. Support for this work has been provided by the lAC (grant 
310394), and the Education and Science Ministry of Spain (grants AYA2007-3E3506, and 
AYA2010-16717). 
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Table 1: Coordinates, radial velocities RV ([km s ^]), projected rotational velocities vsini 
([km s~^]), and atmospheric parameters for the target stars. 



ID 


a(J2000) 


5(J2000) 


RV 


vsini Teff logp ^t" 


Teff logff [A/H] 


6 


logs 










isochrones 


adopted 




Ha 



166 


18;36;07.73 


-23:55:43.4 


-138.8 


15 


7970 


3.10 


2.56 


7700 


3.95 


-1.76 


1.50 


3.30 


2212 


18;36;30.36 


-23:57:12.7 


-153.0 


31 


7643 


3.02 


2.65 


7550 


3.02 


-1.55 


2.67 


2.69 


2297 


18:36:31.12 


-23:51:45.1 


-142.3 


17 


8317 


3.17 


2.48 


7850 


3.17 


-1.83 


2.50 


3.02 


2346 


18:36:31.66 


-23:49:30.2 


-144.5 


17 


7445 


2.99 


2.70 


7400 


2.99 


-1.63 


2.40 


2.81 


2548 


18:36:36.96 


-23:54:33.1 


-142.4 


16 


7366 


2.97 


2.72 


7330 


2.97 


-1.70 


2.10 


2.71 


2562 


18:36:37.68 


-23:57:24.3 


-136.1 


26 


8317 


3.20 


2.48 


8317 


3.20 


-1.61 


2.10 


3.36 


90 


18:36:02.44 


-23:53:38.2 


-148.2 


12 


8227 


3.19 


2.50 


8000 


3.19 


-1.67 


2.10 


3.13 



"From the Pace at al. (2006) empirical relation, by assuming temperatures from isochrones. 



This preprint was prepared with the A AS IM^jX macros v5.2. 



-17- 



Table 2: Derived chemical abundances for the M22 HB stars, and mean RGB stars abun- 
dances from Marino et al. (2009, 2011) in the entire, the s-poor and the s-rich samples. 



ID 


[Fe/Hli^ 


[Fe/H]i 


[Fe/H]„ 


[Fe/H]„ [Na/Fe] 


[Na/Fe] 


[Mg/Fe] 


[Ca/Fe] 


[Ti/Fe]„ 


[Cr/Fe] 


[Ba/Fe] 




LTE 


NLTE 


LTE 


NLTE LTE 


NLTE 


LTE 


LTE 


LTE 


LTE 


LTE 










HB stars 












166 


-1.89 


-1.74 


-1.75 


-1.76 -0.11 


-0.39 


+0.22 


+0.46 


+0.52 


+0.14 


+0.14 


2212 


-1.79 


-1.51 


-1.55 


-1.51 +0.15 


-0.32 


+0.28 


+0.47 


+0.37 


-0.06 


-0.20 


2297 


-2.05 


-1.78 


-1.83 


-1.80 -0.07 


-0.43 


+0.34 




+0.36 






2346 


-1.87 


-1.59 


-1.61 


-1.58 +0.24 


-0.26 


+0.30 


+0.46 


+0.35 


-0.13 


-0.40 


2548 


-1.94 


-1.67 


-1.70 


-1.68 +0.12 


-0.35 


+0.39 


+0.39 


+0.30 


-0.23 


-0.35 


2562 






-1.61 


-1.58 +0.10 


-0.31 


+0.42 




+0.39 






90 


-1.86 


-1.65 


-1.67 


-1.65 -0.08 


-0.39 


+0.34 




+0.32 




< -0.05 


avg. 


-1.90 


-1.64 


-1.66 


-1.63 +0.08 


-0.34 


+0.35 


+0.44 


+0.35 


-0.14 


-0.32 


± 


0.05 


0.0-5 


0.04 


0.05 0.06 


0.03 


0.02 


0.03 


0.01 


0.06 


0.07 




0.10 


0.10 


0.10 


0.10 O.l:-! 


O.(K) 


0.05 


0.01 


o.o:-! 


0.01) 


0.10 










RGB all stars 












avg. 


-1.76 




-1.76 


+0.26 


+0.17 


+0.39 


+0.30 


+0.32 


-0.13 


+0.09 


± 


0.02 




0.02 


0.05 


0.04 


0.02 


0.01 


0.01 


0.02 


0.04 


O'obs 


0.10 




0.10 


0.27 


0.25 


0.11 


0.07 


0.07 


0.08 


0.22 










RGB s-poor (N 


a-poor+Na-rich) 










avg. 


-1.82 




-1.82 


+0.15 


+0.08 


+0.38 


+0.30 


+0.32 


-0.16 


-0.05 


± 


0.02 




0.02 


0.06 


0.05 


0.03 


0.01 


0.01 


0.02 


0.03 


O'obs 


0.07 




0.07 


0.26 


0.24 


0.12 


0.04 


0.06 


0.05 


0.12 










RGB s-rich (N 


a-poor+Na 


-rich) 










avg. 


-1.67 




-1.67 


+0.42 


+0.31 


+0.39 


+0.36 


+0.32 


-0.09 


+0.31 


± 


0.01 




0.01 


0.06 


0.05 


0.03 


0.02 


0.02 


0.04 


0.04 


O'obs 


0.05 




0.05 


0.20 


0.19 


0.11 


0.06 


0.06 


0.08 


0.13 



"The mean values for HB stars have been calculated excluding the star #166. 
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Table 3: Sensitivity of derived abundances to the atmospheric parameters and EWs. We 
report the variations of chemical abundances due to the atmospheric parameters, errors in 
EW measurements (cew)) and the squared sum of the internal errors ((Xtotai)- 





ATeff 


Alog^ 




A[A/H] 


o"ew 


Ctotal 


ATeff 


Alog^ 


A^t 




±170 K 


±0.02 


±0.50 


±0.10 


±8 mA 




±200 K 


±0.15 


±0.30 








internal 






systematic 




[Na/Fe]i 


^0.01 


±0.00 


^0.04 


±0.00 


:f0.08 


0.09 


TO.01 


±0.00 


TO.03 


[Mg/Fe]i 


^0.03 


±0.00 


^0.06 


±0.00 


^0.06 


0.09 


TO.03 


±0.00 


TO.04 


[Ca/Fe]i 


TO.13 


±0.00 


±0.02 


±0.00 


TO. 11 


0.17 


TO.16 


±0.02 


±0.01 


[Ti/Fe]n 


±0.02 


to.oo 


±0.10 


±0.00 


±0.02 


0.10 


±0.03 


±0.00 


±0.06 


[Cr/Fc]i 


^0.13 


±0.00 


±0.02 


±0.00 


TO. 16 


0.21 


TO.16 


±0.02 


±0.01 


[Fc/H]i 


±0.14 


to.oo 


^0.01 


±0.00 


±0.06 


0.15 


±0.17 


TO.02 


to.oo 


[Fc/H]n 


±0.06 


±0.01 


TO.12 


±0.00 


±0.08 


0.16 


±0.08 


±0.04 


TO.08 


[Ba/Fe]„ 


±0.16 


±0.00 


±0.00 


^0.10 


T0.13^ 


0.23 


±0.19 


±0.02 


±0.00 



"Uncertainty introduced by the continuum placement. 
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1 2 3 



B-I 

Fig. 1. — V versus B — I CMD from P. Stetson of M 22 corrected for differential reddening. 
Spectroscopic targets are marked with circles. We superimposed to the CMD the best-fitting 
isochrone from the Padova database (red continuous line), and an isochrone corresponding 
to I.8M0 star (blue dotted line). The inset is a zoom around the HB. 
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Wavelength [k] 



Fig. 2. — Observed spectrum around the Ha region of the star ^^90. The red spectrum is 
the model that best-fits the data (log g'=3.13), while the blue dashed spectra are the models 
with log g higher and lower by ±0.30 dex than the best-fit value. The grey regions are where 
the minimum has been determined. 
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Fig. 3. — Left Panels. [Na/Fe] (top) and [Ba/Fe] (bottom) as a function of Tefr for HB stars 
(circles) and RGB stars (triangles), s-rich and s-poor RGB stars are colored red and blue, 
respectively. The histogram distribution of [Na/Fe] and [Ba/Fe] for HB stars are shown in 
the middle panels, while the red and blue shaded-histograms in the right panels correspond to 
s-rich and s-poor RGB stars, respectively. Measurements for RGB stars come from Marino 
et al. (2011b). Here we have applied a correction of +0.06 dex to the RGB Fe I abundances 
according to the NLTE corrections expected for these stars. 



